Laboratory Search for Spin-dependent Short-range Force from Axion-Like-Particles 

using Optically Polarized 3 He gas 
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The possible existence of short-range forces between unpolarized and polarized spin-i particles 
has attracted the attention of physicists for decades. These forces are predicted in various theories 
and provide a possible new source for parity (P) and time reversal (T) symmetry violation. We use 
an ensemble of polarized 3 He gas in a cell with a 250 fim thickness glass window to search for a force 
from pseudoscalar boson exchange over a sub-millimeter ranges. This interaction would produce a 
NMR frequency shift as an unpolarized mass is moved near and far from the polarized ensemble. 
We report a new upper bound with a factor of 10-30 improvement on the product g 3 gp of the scalar 
couplings to the fermions in the unpolarized mass, and the pseudoscalar coupling of the polarized 
neutron in the 3 He nucleus for force ranges from 10 -4 to 10 -2 m, which corresponds to a mass range 
of 2 x 10~ 3 to 2 x 10 -5 eV for the pseudoscalar boson. This represents the most sensitive search 
that sets a direct limit in the important "axion window" . 
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Possible short-range forces between unpolarized and 
polarized spin-i particles can provide a new source for 
parity (P) and time reversal (T) symmetry violation [l[ . 
Moody and Wilczek Q proposed a force from the ex- 
change of spin-0 bosons which can couple to fermions 
through scalar and pseudoscalar vertices. The scalar 
coupling is spin- independent and depends only on the 
fcrmion density. The pseudoscalar coupling is entirely 
spin-dependent. The resulting spin-dependent short- 
range force (SDSRF) has a Yukawa-type interaction po- 
tential from one boson exchange of the form 
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where r is the unit vector from the unpolarized particle 
to the polarized particle, a is the spin of the polarized 
particle, m p is the polarized particle mass, g s g p is the 
product of couplings of the scalar vertex in the unpolar- 
ized matter and the pseudoscalar vertex of the polarized 
particle, and A is the force range. Such forces may be 
induced by pseudoscalar bosons like the axion [3j , axion- 
like-particles (ALPs) 4] or a very light spin-1 boson [f| 
which are candidates for cold dark matter |6}. Current 
experimental and astrophysical observation restricts the 
axion mass between 1 fj,eV to 1 meV, corresponding to a 
force range between 2 cm to 20 /zm, the so-called "axion 
window" p}. Many ALPs, predicted by stringtheory Q 
and many extensions to the standard model [3, H3] j also 
predict weak forces in this range. 

Several experiments have been performed to search 
for SDSRF using different techniques. Some examples 
include the torsion pendulum [Ill4l3j . neutron bound 
states on a mirror in the Earth's gravitational field [3] , 
and longitudinal and transverse spin relaxation of polar- 



ized neutrons and 3 He A potential of the form 

a ■ f can introduce a shift in the precession frequency 
of polarized particles in the presence of an unpolarized 
mass [TH, similar to that of a magnetic dipole in 
an external magnetic field, oc jl ■ B. First measurement 
using this idea with 3 He [2l[ achieved a sensitivity of 
5 x 10 -3 Hz and restricted the coupling strength close to 
the current limit for force ranges from 10~ 4 to 10~ 2 m 
without any magnetic shielding. In this work, we present 
new results with a factor of 10-30 increase in sensitivity 
which constitute to our knowledge, the most stringent 
laboratory limit on g s g p in the important "axion win- 
dow" . The s pin of the 3 He is dominated by the spin of 
the neutron [22j so this result is directly interpretable 
in terms of the coupling g™. Constraints on T-odd and 
P-odd interactions of the 3 He atom from bounds on the 
electric dipole moment (EDM) of its constituents from 
new physics at high-energy scales are highly suppressed 
due to the Schiff 's screening [23| , and the cancellation of 
the electronic EDM in the ground state of the 3 He atom. 
Our work therefore also represents to our knowledge the 
most sensitive search for T-odd and P-odd interactions 
in the 3 He atom at low energies. 

The apparatus used in this work is based on the de- 
sign of [2l|. We use a 7 amg high-pressure 3 He cell as 
shown in Fig. [TJ which has an optical pumping chamber 
and a target chamber connected by a glass tube. 3 He 
is polarized using spin-exchange optical pumping [2~ij in 
the spherical pumping chamber of radius 4.3 cm. The 
polarized 3 He atoms diffuse into the lower 40-cm long 
cylindrical chamber, which possesses two hemispherical 
glass windows at both ends with a thickness of 250 /im. 

We describe below a number of experimental improve- 
ments compared to [2lj |. which enhance the sensitivity 
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FIG. 1: Schematic of the SDSRF experiment (not to scale). 
The cylindrical polarized 3 He cell is located in a uniform mag- 
netic field. Correction coils (dashed-loop curves) compensate 
for residual holding field gradients. The direction of the laser 
and the holding field are along z. 



to SDSRF. A pair of correction coils is applied to the 
Helmholtz coils to improve the uniformity of the holding 
field. Two identical pickup coils A and B of 2.5 cm di- 
ameter are located next to each other at the same end 
of the 3 He cell. Pickup coil A is mounted below the 
window to measure the precession frequency shift of the 
polarized 3 He nuclei due to SDSRF from the unpolarized 
mass. Pickup coil B is positioned to be insensitive to SD- 
SRF; its signal is used to monitor the holding field drift. 
We subtract the frequencies measured in both coils and 
form f' A = /a — /b for each measurement. The 3 He cell 
position relative to the Helmholtz coils is adjusted to op- 
timize the transverse spin relaxation time measured 
from coil A and B. 

The holding field is tuned to produce a 3 He Larmor 
frequency near 23.8 kHz. We apply a 24 kHz RF pulse 
to tip the spins by a small angle with negligible polar- 
ization loss. The polarized 3 He nuclei induce EMFs in 
the pickup coils which are digitized and recorded. The 
precession frequency is determined first by applying a 
Fourier transform to a signal s(t) in the time domain 
and obtaining the real and the imaginary parts of the 
signal in the frequency domain as R(f) and /(/), respec- 
tively. The Fourier transform is numerically calculated 
using Richardson extrapolation (25| . The total amplitude 
is S(f) = \/R 2 (f) + I 2 {f)- The reference frequency / is 
then varied with a 10 -6 Hz step to locate the maximum 
of S(f), which is the precession frequency [26j |. 

Two samples are used as the unpolarized masses: a 
Macor ceramic mass block of dimensions 34 x 52 x 38 
mm 3 used in [2l| and a liquid mixture of 1.02% MnCl2 
in pure water. These samples are chosen for their dif- 



ferent nucleon densities, low magnetic impurities and 
magnetic susceptibilities, and minimal influence on the 
nuclear magnetic resonance (NMR) measurement proce- 
dure. The paramagnetic salt is added in order to com- 
pensate for the diamagnetism of the water. The magnetic 
susceptibility of this mixture is measured to be < 5% of 
that of pure water. A stepping motor is used to move the 
ceramic mass a distance of 5 cm to 10 /xm from the target 
chamber window. The salt water is stored in a cylindri- 
cal PTFE tank of a radius 34 mm and a length 37 mm. 
Its one end is sealed with a 25 /im flexible PTFE film, 
which contacts the target chamber window. The other 
end of the liquid tank is connected to a flexible PTFE 
tube, which moves the liquid in and out of the tank us- 
ing a nonmagnetic air cylinder actuated by a magneti- 
cally shielded switch. 

In this work, we also improved the analysis method, 
which is described below. We define a mass-in state 
(in) or mass-out state (out) with the mass close to or 
away from the chamber window. Each measurement cy- 
cle employs two states of the mass position in the se- 
quence (in, out, in, out) with a 60 second pause in 
the middle. We apply the analysis algorithm presented 
in [27j to derive the frequency difference between the 
two states and remove any possible bias from linear 
or quadratic time-dependent frequency drifts. Assum- 
ing linear and quadratic time-dependent frequency drifts 
f(t) oc at + bt 2 ±c with a and b being arbitrary constants, 
and +c(— c) the frequency shift depending on the in(out) 
state, the frequency difference between two successive cy- 
cles 1 and 2 is given by 

A/ = —[fin,l ~ 3f out .l + ifin,2 ~ /out, 2] 

= -[(aSt + bSt 2 + c) - 3(a(2St) + b(2St) 2 - c) 

+ 3(a(3St) + b(3St) 2 + c) - (a(A5t) + b(ASt) 2 - c)] 
= 2c (2) 

where St is the measurement time step (the time at the 
beginning of the first step of cycle 1 is taken as zero) and 
/in/ out, 1/2 i s the frequency measured in the pickup coil A 
minus the pickup coil B for cycles 1 and 2, respectively. 
Higher-order algorithms produced the same results. 

The mass in-mass out frequency difference can be mea- 
sured in four different configurations of the apparatus 
corresponding to the directions of the main holding field 
and of the 3 He polarization [2l|, each of which should 
possess the same magnitude of a frequency shift in the 
presence of a nonzero SDSRF proportional to the nu- 
cleon density of the mass. However, for our apparatus, 
two of these configurations possess residual field gradients 
in the sample large enough to lower the spin relaxation 
time T2 and produce complicated line shapes whose fre- 
quency shifts, determined by a peak-finding algorithm of 
the type used in our analysis, are too sensitive to possi- 
ble magnetic systematic effects induced by the mass. We 
therefore consider only two of the four configurations of 
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Sample 


+ -[10" 5 Hz] 


- ±[10~ 5 Hz] 


A/ srf [10- 5 Hz] 


ceramic 


0.6 ±1.3 


-4.6 ±3.1 


2.6 ± 1.7 


salt water 


-3.3 ±0.8 


-1.7 ±5.2 


-0.8 ±2.6 



TABLE I: Data of each configuration of two samples. A/s,p 
is the average frequency difference between mass-in and out 
states and the frequency measured in the pickup coil A minus 
the pickup coil B of each configuration. 



the apparatus which lead to longer spin relaxation time 
T 2 . 

We define the precession frequency in different configu- 
rations of the holding field and the polarization direction 
as Ib,p where the holding field B = ± and the polariza- 
tion direction P — ±. The precession induced by SDSRF 
does not change after reversing the holding field, assum- 
ing the holding field rotation does not change the polar- 
ization of the 3 He. But its precession direction now be- 
comes opposite to the magnetic field-induced precession. 
Any systematic effects depending only on the polariza- 
tion direction do not change after reversing the holding 
field. For different configuration, we can write 



f+± = Jb ± A/p ± A/sDSRF, 
f-T = /s ± A/p - A/ S DSRF 



(3) 



where /p is the magnetic field-dependent precession fre- 
quency including the holding field and the possible effect 
from the magnetic susceptibility of the mass, A/p is the 
polarization-dependent frequency shift and A/sdsrf rep- 
resents the frequency shift due to SDSRF. In this work, 
the two configurations (H — ) and (— ±) which both have 
clean line shapes are considered in determining SDSRF 



A/ SDS RF = i(A/ + _-A/_ 



(4) 



where A/_| and A/ |_ are the frequency differences be- 
tween the mass-in and mass-out states and the frequency 
measured in the pickup coil A minus the pickup coil B of 
each configuration. 

We take 1000 cycles continuously for each configura- 
tion of each sample. The uncertainty in the measured 

t 2 ,. Fig. [2] and 



frequency shift is given by \\J o\_ + 

Table U show the data of two samples and the average 
frequency shift due to SDSRF. The result shows that the 
average frequency difference of the salt water is consistent 
with zero. 

Based on these results, we can constrain the force range 
and the coupling strength. The precession frequency shift 
due to SDSRF for each polarized 3 He nucleus in the tar- 
get chamber can be calculated by numerically integrating 
Eq.([T]) over the unpolarized mass as 



2N f 

Af(z,X,g s9 ;) =— \ 
2irh J vol 



V(r~z)dr 3 



(5) 



where N is the particle number density of the unpolarized 
mass, vol is the total volume of the unpolarized mass 
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FIG. 2: Frequency difference of Macor ceramic(top) and salt 
water(bottom) for 1000 cycles. The average A/ Sr f is obtained 
by fitting the distribution as a Gaussian function. Small boxes 
show the A/ sr f verse cycle. 



and z is the distance from the surface of the mass to the 
polarized 3 He nucleus. The precession signal measured 
by the pickup coil is s(t) oc J (B CO ii ■ dM/dt)dr 3 where M 
is the magnetization vector of 3 He and B co u is the field 
profile of the pickup coil which can be derived by using 
the reciprocity theorem [28]. The signal is written as 

s{t) = C{j -B x ■ (Jo + A/) sin(27r(/ + Af)t) 

Jd 

+ B y ■ (f + A/) cos (2tt(/ + Af)t)dz} (6) 

where /o = "fBo/2ir is the Larmor frequency, 7/277 = 
—3.24 Hz/mG is the gyromagnetic ratio of 3 He, d is the 
thickness of the cell window and C is a constant. Ap- 
plying the Fourier transformation, the power spectrum 
of the signal can be calculated as, 

P(f') = C{( / (/o + Af)BJ(f ± A/ - f')dz) 2 



+( / (/o + Af)B y S(f ± A/ - f')dz) 2 }. (7) 



(8) 



The average frequency observed by the pickup coil is 
If'P(f')df 



JP(f')df 



The frequency shift due to SDSRF is Af(\,g s g%) = 
f — fa- Using the measured frequency difference in Ta- 
ble HI the constraint on the coupling strength and the 
force range is determined as shown in Fig. [3] where the 
dark gray area was ruled out by previous measurements. 
The dotted curve is from [Ti| and the dash-dotted curve 
is from [l7[ ■ The dashed (solid) curve is the constraint of 
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FIG. 3: Constraints on the coupling strength g s gp as a func- 
tion of the force range A and the equivalent mass of the ALPs. 
The dark gray area is the region excluded by previous works. 
The dotted curve is from [10] and the dash-dotted curve is 
from [13]. The dashed (solid) curve is the constraint of the 
salt water (ceramic) sample within one standard deviation. 

the salt water (ceramic sample) within one standard devi- 
ation. The measured frequency difference of the ceramic 
sample due to SDSRF is consistent with zero within 1.5 



standard deviations. Our new results improve the con- 
straint on SDSRF from the current limit in the range of 
10~ 4 to 10~ 2 m by a factor of 10-30, which corresponds 
to a mass range of 2 x 10~ 3 to 2 x 10~ 5 eV for the pseu- 
doscalar boson involved. This work represents the most 
sensitive search that sets a direct limit in the important 
"axion window" . 

Several methods can be employed in the future to fur- 
ther improve the sensitivity using polarized 3 He. Obvious 
paths for improvement of the measurement include new 
magnetic holding field systems with better field unifor- 
mity and magnetic shielding, a smaller 3 He cell with a 
lower pressure and thinner windows, unpolarized mass 
samples with higher fermion densities and lower mag- 
netic susceptibilities, and a 129 Xe comagnetometer. With 
these changes, we conclude that a factor of 10-100 im- 
provement in the constraints of the coupling strength in 
the force range of 10~ 4 to 10~ 2 m is possible. 
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